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We study the mechanism of dynamical nuclear spin polarization by hyperfine interaction in spin- 
blocked double quantum dot system. We calculate the hyperfine transition rates and solve the master 
equations for the nuclear spins. Specifically, we incorporate the effects of the nuclear quadrupole 
coupling due to the doping-induced local lattice distortion and strain. Our results show that nuclear 
quadrupole coupling induced by the 5% indium substitution can be used to explain the recent 
experimental observation of missing arsenic NMR signal in the spin-blocked double dots. 
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In many metals and semiconductors, nuclei and con- 
duction electrons form a coupled nonlinear spin system 
through the hyperfine interaction. Such a coupling is well 
studied in bulk experimental phenomena such as dynam- 
ical nuclear spin polarization Both electron and nu- 
clear spins in semiconductor heterostructures have been 
proposed as candidates for qubits in a scalable quantum 
information processor. There are also a variety of spin- 
tronic devices where electron spins provide novel func- 
tionalities. For both quantum computation and spintron- 
ics, understanding the electron-nuclear spin interaction 
and the resulting coupled dynamics is crucial because 
of its important roles in spin manipulation and decoher- 
ence. Indeed, in the past decade there have been ex- 
tensive experimental and theoretical studies of electron- 
nuclei h3rnerfine interaction in semiconductor nanostruc- 
tures In this Letter, we show that nuclear 

quadrupole coupling can be a critical factor in this cou- 
pled dynamics, thus needs to be carefully accounted for 
in the study of spin-based quantum computers and spin- 
tronic devices. 

Dynamical nuclear spin polarization has recently been 
demonstrated in a spin blocked semiconductor double 
quantum dot made of Gao.gsIno.osAs Coulomb 
interaction and Pauli principle dictates that in a double 
quantum dot two-electron singlet and triplet states are 
split, so that proper voltage offset and bias between the 
double dot leads to significant suppression in the tun- 
nel current due to occupied triplet state (thus the so- 
called spin blockade regime, which has been suggested 
for single spin detection, a crucial component of quan- 
tum information processing) 0- One way to lift this 
blockade is to apply an appropriate magnetic field, so 
that electron-nuclear spin hyperfine coupling can facil- 
itate electron spin-flip transitions, which in turn leads 
to dynamical nuclear spin polarization as was shown in 
Ref. 0. The reported experimental observations reveal 
several quite unexpected phenomena [^, one of which 
is that while nuclear magnetic resonance (NMR) signals 
were observed at the ^^a and ^^Ga frequencies, no com- 
parable response has been found from the ^^As nuclei. It 
seems that nuclear spin polarization cannot be built up in 



^^As nuclei even though they have the largest concentra- 
tion and the largest hyperfine coupling with the electrons 
among the three nuclear species. 

In this Letter we study the mechanism of dynamical 
nuclear spin polarization by hyperfine interaction in a 
spin-blocked double quantum dot, taking into account 
electron tunneling and nuclear quadrupole coupling. We 
show that the absence of ^^As nuclear spin signal is due to 
the local strain created by the presence of Indium atoms. 
Specifically, the substitution of the 5% In for Ga breaks 
the crystal symmetry and induces a local lattice distor- 
tion [3, 01 . The charge redistribution generates static 
electric field gradients (EFG) at the surrounding lattice 
sites, which lead to finite quadrupole coupling for the 
nuclei 0| . The randomly distributed In atom only sub- 
stitute one of the nearest neighbors of As atoms, while 
the nearest neighbors of Ga atoms are always four As 
atoms. Therefore the local electric field gradients at the 
Ga nuclei sites are on average much smaller than at the 
As sites. We show that this difference of EFG at the 
locations of As and Ga nuclei is the major origin of the 
unexpected absence of As nuclear spin signal. 

We first note that the static quadrupole interaction 
alone cannot generate any transition among the spin 
eigenstates because there is no energy relaxation mecha- 
nism. This is different from the spin-lattice effect, where 
the EFG is generated by the lattice oscillations. There 
quadrupole spin relaxation can be realized by absorb- 
ing and emitting a phonon simultaneously (i.e., a Raman 
process) |l2j . In the system we consider, the nuclei cou- 
ple with the electrons, which in turn couple to the leads 
(reservoir) by tunneling, so that energy conservation can 
be satisfied through charge exchange between the dot 
and the leads. In the following we first find the nuclear 
spin eigenstates including the static quadrupole interac- 
tion, then evaluate the first-order transition rates among 
these eigenstates due to hyperfine coupling. 

The static Zeeman and quadrupole Hamiltonian for a 
single nuclear spin is 
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where lo is the nuclear spin Lamour frequency and v is 
quadrupole couphng strength that is proportional to the 
EFG at the locations of the nuclei, z is the direction of 
external magnetic field and Z is the principle axis of the 
largest EFG We have neglected the asymmetrical 

effect of quadrupole interaction. It was determined ex- 
perimentally that j^As ~ 2.5 MHz and v an i s much less 
than i/As in Ina;Gai_a;As for x up to 0.01 For a typ- 
ical external field i?ext (0.8 Tesla) applied in the experi- 
ment 0, we can assume i/ca wca ~ 10.4 and 8.2 MHz 
for ^^Ga and ^^Ga nuclei. Since our intention is to study 
the relation of As and Ga nuclear polarization, wc will 
not differentiate ^^Ga and ^^Ga further in this Letter. 
The small I'ca allows eigenstates (i=-3/2,-l/2,l/2 or 
3/2) of Hn for Ga nuclei to be expressed perturbatively 
in terms of the Zeeman spin eigenstates \m). 
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where h = -v/Si^GaSinflcos^/wGa and c = 
-\/3j^GaSin^6'/4a;Ga- d is the angle between z and Z 
axes. In these calculations x axis is chosen to be in 
the plane of z and Z. These new states \(f)i) are not 
eigenvectors of Iz anymore, though their average spin 
polarizations {4ii\iz\4'i) =i + Oiiy^^/uj'l.J. 

Transitions can be induced by the hyperfine interaction 
among the new energy eigenstates of the Ga nuclei. The 
hyperfine Hamiltonian is 
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Here iV ^ 10^ is the number of nuclei in a quantum 
dot; A is the hyperfine coupling constant; ipCRi) is the 
electron envelope wavefunction at the ith nucleus. Using 
the electron density at nuclear sites 0] we estimate A « 
46 ^eV and 40 /ieV for As and Ga nuclei respectively. 
Hamiltonian Q describes the interaction of an electron 
with an ensemble of N nuclei in the dot. 

To calculate the transition rates given in Fig. 1 we must 
study how the system couples to the environment to sat- 
isfy energy conservation. In the present case electron 
exchange between a dot and the leads is the most impor- 
tant energy relaxation mechanism [isj . For triplet state 
I ii) the electron in the 1st (2nd) dot only couples 
to the left (right) lead. The electron tunnel Hamiltonian 
can be expressed as 

He = ^efcc+ Cfe^+eo4rfi+^ [tfeC+ rf^ + H.C.] , (4) 

k k 

where eg is the renormalized single particle energy 
level. Since the Hamiltonian is quadratic, ex- 
act solution of the retarded Green's function G^^^* = 
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FIG. 1: Schematic of hyperfine induced transitions among the 
energy eigenstates of Hamiltonian Hn. The downward arrows 
represent simultaneous flip-fiop scattering of nuclei and the 
electron in a triplet state. F and G are the nuclear spin 
relaxation rates induced by the effective hyperfine magnetic 
field of the electron {-^FSz). 



—i9{t){{di{t), d^{0)}) can be obtained by using the equa- 
tion of motion method. The density of state (DOS) fol- 
lows from /9T(e) = — ImGj'j^*(e-|-z0~'')/7r. It is straightfor- 
ward to show that the electron triplet state is broadened 
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where F = 27r |ffep(5(e — e^) is the level broadening 
for the triplet state. For the electron in the singlet state, 
the DOS should be modified by a weak inter-dot cou- 
pling. However, the major energy relaxation channel is 
still the lead-dot tunnel coupling, so we can approximate 
the sing let DOS with Eq. © as weU. 

The first term in Eq. (jS)) does not involve electron spin 
flip. However, it can induce transitions among the energy 
eigenstates due to mixture of different spin eigenstates in 
Eq. The transition rate between state \(f>i) and {(pj) 
due to this term is given by the Fermi golden rule, 
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where Ei and Ej are the eigenenergies of i/„. The ab- 
sorption power P is 



P(A)= ^p(6,)/^(6/)<5(e/-e, + A), 
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where p(e) is the DOS of either the electron singlet or 
the triplet state. Using Eq. jnj, the transition rates in 
Fig. Ill) can be calculated: F = Qh"^ A^ / {N^fiT) and G 
= 6c^A^/{N^h,r). In deriving these expressions we have 
neglected the nuclear energy difference of the initial and 
final state since they are much smaller than the width of 
electronic energy broadening. 

The calculation of the transition rates due to the hy- 
perfine flip-flop terms in Eq. 0, which correspond to the 
polarized triplet to singlet electron transitions, is more 
complicated. When the two electron states are degen- 
erate, the hyperfine flip-flop transition is on resonance. 
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leading to nuclear spin polarization. However, the tran- 
sition should quickly becomes off resonance due to the 
effective nuclear magnetic field Bjq experienced by the 
electrons. The effective nuclear magnetic fields can be 
estimated as 
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where 7as « -2.8 Tesla and 7Ga « -1.3 Tesla As 
and Ga nuclei have independent spin temperatures be- 
cause their mutual flip-flop process is largely suppressed 
due to different magnetic moments. Equilibrium in nu- 
clear spin polarization will thus be established in each 
of the nuclear species independently. Bj^ can be quite 
large for even moderate nuclear polarization. The en- 
ergy difference of the singlet state and triplet state due 
to can be approximated by Ast = — .9*M-B-^n as- 
suming the energy of the singlet state is independent of 
magnetic field. The flip-flop transition rates of hypcrfinc 
interaction from state to state is thus 
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where /t is a multiplication factor due to cotunneling and 
the probability of the triplet state occupation (electrons 
can tunnel into three triplet states and one singlet state, 
but only one triplet is relevant to nuclear polarization). 
/t should be multiplied by another factor 0.5 because 
there are two dots and the probability of nuclear spin- 
flip process in each of the dots is 50%. We estimate 
/t w 0.05. We then flnd that D in Fig. is given as 
D = fTA^T/N^h{AlT + r2). 

The quadrupole interaction of As nuclei is much 
stronger than that of Ga [z^as ~ 2.5 MHz is in the same 
order as the Zeeman energy ujas ~ 5.9 MHz (at 0.8 Tesla)] 
and leads to a complete mixing of the spin eigenstates, 
so that the perturbative wavefunctions [Eq. JSJ] for Ga 
nuclei have to be replaced with numerical solutions: 
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The expectation value of Iz of these state is {4>i\iz\(t)i) = 
X^m [""imp"*- For example, we find that when a = 
VAs/^As = 0.3, {(p^\iM^) = 1-419, 0.558, -0.528 and - 
1.450 for i = 3/2, 1/2, -1/2 and -3/2 respectively All 
the matrix elements of Iz and /_ can be easily obtained 
from Eq. Hl()|l , and the transition rates can be calculated 
numerically using Eq. and ©. One should keep in 
mind that the illustration of the various transitions in 
Fig. is only applicable to the perturbative analysis for 
Ga nuclei, while for As nuclei all transitions are possible. 
However the concept of polarization {D) and depolariza- 
tion (F and G) processes illustrated in Fig. should 
still be valid. 




FIG. 2: Ratio of Ga and As nuclei polarization as a function 
of the relative quadrupole interaction strength a {— voa./i'As) 
for r = 10 /ieV (panel a) and F = 30 /xeV (panel b). F is the 
width of level broadening due to electron tunneling between 
the lead and the dot. The three curves in each graph represent 
three different As quadrupole coupling strengths vab/uiab = 
0.1, 0.3 and 0.5. 6 has chosen to be 45 degrees for both As 
and Ga nuclei. 



After obtaining all the transition rates W and W 
we can proceed to construct the master equations [l^ 
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nuclear spins in state In Eq. Hll|l we have neglected 
the spin-lattice relaxation, which is very weak at liquid 
Helium temperature. Here we can use the spin temper- 
ature approximation because the nuclear spin build-up 
time (« 100 s) is much longer than the transverse spin 
relaxation T2 (« 10"'' s). We solve for the steady-state 
polarization: dpi/dt = 0. These equations are highly 
nonlinear, as the equilibrium polarizations depend on the 
transition rates, while to calculate the transition rates 
one must know the effective nuclear magnetic field, 
which itself depends on the nuclear spin polarization. We 
solve the set of equations self-consistently. The average 
nuclear polarization is P = '}2iPi{'^i\^A4'i) ~ {Iz)/I- We 
can then calculate the effective nuclear magnetic fields 
Sn with Eq. ©. 

Figure|21shows the calculated ratio of Ga and As nuclei 
polarization as a function of a, which is taken as a free pa- 
rameter in our calculations. We consider As quadrupole 
effects with three different strengths. i^as/^^As = 0.3 cor- 
responds to j'As ~ 1-8 MHz. In this case the Ga nuclear 
polarization is two orders of magnitude greater than that 
of As nuclei when a « 0.01. This difference can quali- 
tatively explain the experimental observation that the 
As NMR signal is missing in measuring the current os- 
cillation period as a function of the AC magnetic field 
frequency. Wc also find the ratio is very sensitive to elec- 
tron level width F. Figure [3 shows that the polarization 
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FIG. 3: Ga nuclear spin polarization as a function of the 
relative quadrupole interaction strength a for F = 10 fieV 
(panel a) and F = 30 fieV (panel b). All the parameters are 
the same as those used in Fig. 2. 



to static quadrupole interaction induced by lattice distor- 
tions. Our results can thus be used to explain the recent 
NMR experiments conducted in such system. Our cal- 
culation suggests both caution and promise. Doping is 
widely used in semiconductors to tailor electronic proper- 
ties. However, as we point out in this Letter, there can be 
unexpected side effects in material properties, such as in- 
creased nuclear spin relaxation in the present case. Con- 
versely, controlled doping can also be used to differentiate 
parts of a system (such as Ga and As nuclei here), so that 
selective operations become possible. How to utilize the 
additional control provided by doping while overcoming 
its negative effects can be critical to future quantum in- 
formation processors and/or spintronic devices. 

We acknowledge helpful conversations with K. Ono, S. 
Tarucha, and S. Das Sarma. We thank partial financial 
support by ARDA and ARO. 



ratio in the case of F = 30 jueV is roughly i of that in the 
case of F = 10 /xeV. Since D cx F/(A|j, F^), the triplet 
to singlet transition rate decreases much faster when F is 
smaller, resulting in larger polarization ratio. 

Figure |31 shows how Ga nuclear spin polarization 
changes with the relative quadrupole coupling strength 
a. It is obvious from the two panels that Pqh increases 
rapidly as a decreases from 0.1 to 0.01. The Ga polariza- 
tion saturates to unity for smaller a, and larger F leads 
to larger Ga polarization. 

The physical picture of the master equation calcula- 
tion is actually quite straightforward. Both Ga and As 
have a polarization rate and a depolarization rate. The 
As nuclear spin depolarization rate is much larger than 
that of Ga due to the quadrupole coupling. Thus there 
exists a regime of polarization rate where nuclear spin 
polarization can be built up in Ga but not in As. 

In our study As depolarization rates are greatly en- 
hanced because the energy levels of the dots arc broad- 
ened through tunnel coupling to the leads, and that the 
singlet-triplet transition becomes off resonance once the 
polarization of Ga nuclei is established. In previous ex- 
periments concerning dynamical nuclear polarization, op- 
tical pumping creates a constant polarization rate for the 
nuclei with the help of hyperfine interaction; and the re- 
laxation process due to carrier recombination overshad- 
ows the effect of quadrupole interaction, so that As po- 
larization can still be built up . 

In conclusion, we have studied the dynamical nuclear 
polarization mechanism in spin-blockade double-dot sys- 
tem. Specifically, we calculate the hyperfine induced po- 
larization rates and depolarization rates among the eigen- 
states of the nuclear spins. We show that the average spin 
polarization in steady state of different nuclear species 
(Ga and As) could differ by two orders of magnitude due 
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